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^ ■ Abstract 

o. 

"^ , QCD coherence effects are studied based on measurements of correlations of parti- 

D I cles with either restricted transverse momenta, pT < p^^^ , where pT is defined with 

p ■ respect to the thrust axis, or restricted absolute momenta, p = |p| < p™*, using 

O , about four million hadronic Z decays recorded at LEP with the OPAL detector. The 

correlations are analyzed in terms of normalized factorial and cumulant moments. 
The analysis is inspired by analytical QCD calculations which, in conjunction with 
S^ . Local Parton-Hadron Duality (LPHD), predict that, due to colour coherence, the 

J^ I multiplicity distribution of particles with restricted transverse momenta should be- 

come Poissonian as p^^^ decreases. The expected correlation pattern is indeed 
observed down to p^^ ~ 1 GeV but not at lower transverse momenta. Further- 
more, for p™* ^0 GeV a strong rise is observed in the data, in disagreement with 
theoretical expectation. The Monte Carlo models reproduce well the measurements 
at large p™* and p*^"* but underestimate their magnitudes at the lowest momenta. 
The e+e~ data are also compared to the measurements in deep-inelastic e"'"p colli- 
sions. Our study indicates difficulties with the LPHD hypothesis when applied to 
many-particle inclusive observables of soft hadrons. 
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1 Introduction 

At high energies, the annihilation process e"'"e~— > hadrons proceeds through the creation 
of a highly virtual primary quark and anti-quark which initiate a cascade of partons 
through successive parton emissions. The evolution of such a parton cascade is well 
understood in perturbative Quantum Chromodynamics (QCD) for virtualities, Q^, of 
the daughter partons larger than Qg- Here Qo is a virtuality cut-off below which the 
strong coupling constant becomes large and perturbative methods cease to be valid. In 
comparisons to experimental data, Qo "was found to be of the order of hadronic masses 
(Qo ~ few hundred MeV). 

A fundamental property of a QCD cascade, which follows from the non-Abelian struc- 
ture of QCD, is colour coherence. This induces an angular ordering of subsequent emis- 
sions in the branching process [1] which restricts the phase space for each subsequent 
parton in the cascade. Angular ordering has important consequences of which we men- 
tion only a few (see [2-5] for a comprehensive review). Compared to a cascade without 
angular ordering, the single-parton inclusive distribution is suppressed for soft, or low- 
momentum, particles (the "hump-back" plateau) [6], the mean parton multiplicity evolves 
less rapidly with increasing jet energy, and the rapidity distribution becomes flat and en- 
ergy independent for partons with very small transverse momenta [3]. 

It is remarkable that inclusive characteristics of hadrons measured in a variety of hard 
processes indeed show a behaviour similar to that expected from perturbative parton-level 
calculations [3-5]. This indicates that perturbative QCD effects and colour coherence in 
particular leave their imprint on the hadronic final state even for quantities which are 
not infra-red safe, such as particle multiplicity. The hypothesis of Local Parton-Hadron 
Duality (LPHD) [6] embodies these observations. According to LPHD, parton-level QCD 
predictions are applicable to sufficiently inclusive hadronic observables without the need 
for a hadronisation phase: hadronic spectra are proportional to those of partons if the 
cut-off Qq is decreased towards a small value of the order of A, the QCD-scale. 

Within the LPHD picture, perturbative QCD calculations have been carried out in the 
Double Leading Logarithmic Approximation (DLLA) or in the Modified Leading Loga- 
rithmic Approximation (MLLA) which includes terms of order y^os in the strong coupling 
constant [2,3]. The DLLA calculations neglect energy-momentum conservation in gluon 
splittings which is partly taken into account in the MLLA. Although analytical calcula- 
tions provide much valuable physical insight, they have often to be considered as quali- 
tative. More quantitative results are obtained from parton-shower Monte Carlo models, 
the physics implementation of which strongly resembles the analytical calculations, but 
which impose energy-momentum conservation and include the complete parton-splitting 
functions. 

In spite of its success with single-particle inclusive spectra, earlier studies have shown 
that the applicability of LPHD is less evident for the moments of single-particle densities 
at HERA energies [7] and angular correlations at LEP [8] and at HERA [9] . It is therefore 
of considerable importance to further test multiparticle aspects of perturbative QCD 
predictions, which are sensitive to colour coherence, in conjunction with LPHD. 



Sensitive studies of colour coherence were suggested in [10] within DLLA calculations 
and using factorial moment and cumulant techniques. In a QCD cascade the presence 
of one gluon enhances the probability for further gluon emissions, causing positive cor- 
relations. The multiplicity distribution of partons in a jet is therefore generally broader 
than a Poisson distribution (which corresponds to uncorrelated production) and obeys 
asymptotic KNO-scaling [11]. However, it was pointed out in [10] that, due to colour 
coherence, gluons produced with bounded transverse momenta, px < Px"*, where pt is 
defined with respect to the primary parton in a jet, become, for small p^^^, independently 
emitted from the primary parton. This implies that their multiplicity distribution be- 
comes Poissonian, analogous to that of soft photons radiated from a charged particle in 
QED. In contrast, for gluons with bounded absolute momenta, p = |p| < p'^'^^^ for which 
the angular ordering constraint is less important, the distribution remains non-Poissonian 
even for very small p'^"*. 

The DLLA analytical predictions were first tested in deep-inelastic e"'"p scattering 
by the ZEUS experiment at HERA [12] using a sample of ~ 7500 high-Q^ events. Be- 
cause of low statistics, factorial cumulants were not studied. The factorial moments were 
measured in the current region of the Breit frame [13]. From the significant discrepan- 
cies between data and both the DLLA calculations and parton-level Ariadne Monte 
Carlo [14] expectations, the authors conclude that the LPHD hypothesis is strongly vi- 
olated for many-particle observables. Monte Carlo models, which include hadronisation 
effects, reproduce the correlation pattern of the hadronic final state, although sizeable 
discrepancies remain for small values of p^"* and p ™* . 

In this paper, we report the first results in e"'"e~ annihilation on factorial moments 
and cumulants for hadrons with restricted transverse and absolute momenta in a jet. 
The measurements are based on a data sample of about four million Z hadronic decays 
recorded with the OPAL detector at the LEP e"^e~ collider at CERN. 

2 Analysis 

The calculations in [10] use factorial moments and cumulants known to provide a sensitive 
tool to probe multiparticle correlations [3,4,15,16]. 

The normalised factorial moment of order g in a region of phase space of size Q is 
defined as 

F,(^) = (n(n-l)••■(n-g+l))/(n)^ q ^ 1. (1) 

Here n is the number of particles in Q and the angle brackets (■ ■ ■) denote the average 
over events. For uncorrelated particle production within Q one has Fg = 1 for all q. The 
factorial moments describe many-particle distributions via the relation {n{n — 1) ■ ■ ■ {n — 
g + 1)) = J^pq{pi, . . . ,pq) riLi'^P* between the unnormahsed factorial moments in a 
region Q and the inclusive g-particle densities, Pq{pi, ■ ■ ■ ,Pq), of particles with momenta 

Pi- 

The normalised factorial cumulants, Kq{Q), or cumulants for short [17-19], are related 



to the factorial moments Fq{Q) through the following relations: 

K2 = F2-1, K3 = F3-3F2 + 2, K^ = F^ - AF;, - 3 F^ + 12 F^ - 6. (2) 

By construction, Kg is a measure of genuine multiparticle correlations: Kg, representing 
the correlation function averaged over the region Q, vanishes whenever any one of the q 
particles is statistically independent of the others. For uncorrelated particle production, 
or Poissonian emission, within Q one has Kg = for g > 1. 

The normalised factorial moments of the multiplicity distribution of gluons which are 
restricted in either transverse momentum px < p^^^ (cylindrically-cut phase space) or 
absolute momentum p < p'^^^ (spherically-cut phase space) are predicted to have the 
following qualitative behaviour [10]: 

F rr,™M ~ 1 ^ ^(^~^^ ln(p^"VQo) 

and 

Fg(p™') ~ Ci{q) > 1 for p™* ^ GeV, (4) 

where E is the energy of the initial parton, and transverse momentum is defined with 
respect to its direction. Here and below the C-functions are g-dependent constants. The 
cumulants (0) are predicted to behave as 

*'<^^"*' - (tot)'" "" ^"^«°- ''' 

Equation (j3)) shows that the lower-order cumulants dominate at small p™*. The spherically- 
cut cumulants are predicted to behave similarly to factorial moments, Eq. (^: 

ir,(p™*) ~ C2(g) > for p^"' -^ GeV. (6) 

Equations (jSl)-® illustrate the different influence angular ordering has on the multiplic- 
ity moments and correlations. Cylindrically-cut moments show positive correlations but 
approach the Poisson limit as p™^ approaches Qo- On the other hand, for soft gluons 
with limited absolute momenta, p < p^^^^ the multiplicity distribution remains broader 
than a Poisson distribution for any (small) value of p'^"*. 

In [10], the analytical results have been tested at parton level using the shower Monte 
Carlo program Ariadne [14]. The moments Fq{p^^^) indeed show the expected decrease 
for small values of the cut, p^"* < 4 GeV; however, they do not fully reach the Poisson 
value for p™* — > Qq but saturate at values somewhat larger than one. On the other hand, 
Fg{p^^^) moments attain a maximum around p'^"* ^ 2 GeV, and then decrease towards 
finite values much above unity for p'^"* -^ GeV showing that in spherically-cut phase 
space there is no Poisson regime. 



3 Experimental details 

3.1 The OPAL detector 

The OPAL detector, operated from 1989 to 2000 at LEP, is described in detail else- 
where [20]. The results presented here are mainly based on the information from the 
tracking system, which consisted of a silicon microvertex detector, an inner vertex cham- 
ber, a jet chamber with 24 sectors each containing 159 axial anode wires, and outer 
^-chambers to improve the z coordinate resolution^. The tracking system was located 
in a 0.435 T axial magnetic field and measured p_L, the track momentum transverse to 
the beam axis, with a precision of ((Xp^/px) = a/(0.02)2 _|_ (0.0015 p^)^ (p_L in GeV) for 
Icos^l < 0.73. 

3.2 Data and Monte Carlo samples 

This analysis is based on a data sample of approximately 3.9x10^ hadronic Z decays 
collected with the OPAL detector between 1991 and 1995. About 91% of this sample 
was taken close to the peak of the Z; the remaining part has a centre-of-mass energy, ^/s^ 
within ±3 GeV of the Z peak. 

Further selection criteria are based on a hadronic event selection procedure described 
in detail in [21]. For each event tracks were accepted only if they had at least 20 measured 
points in the jet chamber, the first hit closer than 70 cm to the beam axis, the measured 
closest distance to the e^e~ collision point less than 5 cm in the plane perpendicular to the 
beam axis and less than 40 cm along the beam axis, p_L > 0.15 GeV, and | cos^| < 0.94. 

The event was then required to have at least five tracks, a momentum imbalance, de- 
fined as the magnitude of the vector sum of momenta of all charged particles, below 0.4y^, 
a total energy of the tracks (assumed to be pions) greater than 0.2^8, and | cos^thrl < 0.9, 
where 6'thr is the polar angle of the event thrust axis with respect to the beam direction 
calculated using all tracks as well as electromagnetic and hadronic calorimeter clusters. 
These criteria provide rejection against background from non-hadronic Z decays, two- 
photon and beam-wall interactions, beam-gas scattering, and ensure that the event is 
well contained inside the detector. A total of about 2.9 million events remain after the 
selection has been applied and are used for further analysis. 

The kinematic variables used in the analysis are defined with respect to the event 
thrust axis. To remain consistent with the theoretical calculations and with similar 
measurements in deep-inelastic scattering, factorial moments and cumulants are calcu- 
lated for all charged particles in a single event-hemisphere, defined by a plane per- 
pendicular to the event thrust axis. The particles are assigned positive rapidity, y = 
0.5 ln[{E + pi^)/{E — Ph)], with E and Pl the energy (assuming the pion mass) and longi- 

^OPAL uses the right-handed coordinate system defined with the positive z along the direction of the 
e~ beam and the positive x axis pointing towards the centre of the LEP ring, r is the coordinate normal 
to the beam axis, (p the azimuthal angle with respect to the x axis, 9 the polar angle with respect to the 
z-axis. 



tudinal momentum component of the particle, and a single randomly chosen hemisphere 
for each event is used in the analysis. 

To correct the measured factorial and cumulant moments for the effects of detector 
response, initial-state radiation, resolution and particle decays, we apply the correction 
procedure adopted in our earlier studies [22,23]. Two samples of more than three mil- 
lion multihadronic events each were used, generated with the Jetset 7.4/Pythia 6.2 
Monte Carlo model [24]. The first sample does not include the effects of initial-state 
radiation, and all particles with lifetimes longer than 3xl0~^° s were considered to be 
stable. The generator-level factorial moments, Fq{Q)gen, are calculated directly from the 
charged particle multiplicity distributions of this sample without any selection criteria. 
The second sample was generated including the effects of finite lifetimes and initial-state 
radiation and was passed through a full simulation of the OPAL detector [25] . The cor- 
responding detector-level moments, -Fg(r2)det; are calculated from this set using the same 
reconstruction and selection algorithms as used for the measured data. The corrected 
moments are then determined by multiplying the measured ones by the correction factors 
Uq{fl) = Fg(r2)gcn/i^g(fi)det- The correction factors vary between about 0.85 and 1.2. 

As systematic uncertainties, we include the following contributions: 

• The statistical error on the correction factors Uq{Q): due to the finite statistics of 
the Monte Carlo samples these are comparable to those of the data. 

• Track and event selection criteria variations as in [22,23]. The moments have been 
computed changing in turn the following selection criteria: the first measured point 
was required to be closer than 40 cm to the beam, the momentum was required to 
be less than 40 GeV, the track polar angle acceptance was changed to | cos 6*1 < 0.7. 
The changes in the corrected moments when the analysis is performed with these 
cuts were taken as systematic uncertainties. These changes modify the results by 
no more than a few percent in the smallest phase space regions and do not affect 
the conclusions. 

• Resonance decays: We have repeated our calculations with the Pythia 6.2 Monte 
Carlo model where no decays of resonances were allowed. The difference between 
these calculations and those based on the sample generated including the resonance 
decays are taken as systematic uncertainties and do not exceed 2%. 

• Two-jet selection criteria: For comparison with ZEUS data [12] where the rate of 
hard jet production is lower than in e^e~ annihilation at LEP, we have calculated 
the moments for two-jet events selected via a thrust value cut as given in Sect. 14.31 
Therefore, for the results presented in Sect. 14.31 onlv. we have repeated calculations 
for two-jet events using the Durham jet finder [26]. We apply this algorithm with 
the jet resolution parameter ?/cut = 0.03, shown [27] to result in well separated jets 
while still yielding reasonable event statistics. The changes in the results based on 
the two selection methods are taken as systematic uncertainties. These results agree 
to within 7%. 



• Herwig based correction factors Uq-. The correction factors Uq{Q) were derived from 
samples generated with the Herwig Monte Carlo [28]. The differences compared 
to Pythia are taken as systematic uncertainties and do not exceed 10%. 

The total errors have been calculated by adding the systematic and statistical un- 
certainties in quadrature and are therefore correlated bin-to-bin. It was further verified 
that our conclusions remain unchanged when events taken at energies off the Z peak are 
excluded from the analysis. 

The data are compared to model predictions calculated using the following Monte 
Carlo generators: 

• Pythia version 6.2 [24] with the parton shower followed by string hadronisation. 

• Pythia version 6.2, as above, but including the effect of Bose-Einstein correlations. 
These are simulated using the BE32 algorithm [29] implemented in PYBOEI. In a 
previous OPAL study of higher-order cumulants [23] it was shown that this model 
accounts simultaneously for the magnitude and bin-size dependence of cumulants of 
like-sign as well as of all-charge multiplets in one- to three-dimensional phase space. 
Here we use a Gaussian parametrisation with PYBOEI and QCD/fragmentation pa- 
rameters from [30]. 

• Ariadne version 4.1 [14] with the colour dipole model for the parton shower fol- 
lowed by fragmentation as in Jetset/Pythia. 

• Herwig version 6.3 [28] with a parton shower followed by cluster fragmentation. 

Each Monte Carlo sample consists of more than three million events. The simulation 
parameters of the Jetset/Pythia and Herwig models have been tuned to OPAL data 
in [31]. The parameters of Ariadne and recent changes for the Herwig parameters^ are 
given in [32] . The errors of the Monte Carlo predictions are comparable to those of the 
data. 



4 Results 

4.1 Factorial moments 

Fig. n shows cylindrically cut factorial moments Fq of order g = 2 to 5 as a function 
of Pt"*.^ With decreasing p^"*, the moments decrease towards a minimum at a common 
value of p™* ~ 1 GeV but remain larger than unity, the Poisson value. The observed 
deviation of the pr^"*-dependence from the Poissonian behaviour for large p^"* values agrees 
qualitatively with the theoretical expectation discussed in Sect. |21 However, for smaller 

^Here we used PSPLT(l) = 0.6 and CLMAX== 3.6 GeV instead of 1.0 and 3.35 GeV. 
■^Thc numerical values of the data on factorial moments and cumulants will be made available in the 
Durham HEP Database, http : //durpdg . dur . ac . uk/HEPDATA . 



p^^ values the moments rise strongly, in clear disagreement with the perturbative QCD 
result for partons. Figure ^ suggests that the predicted Poisson limit for soft gluons is 
masked by strong hadronisation effects as p™* — »• Qo- The drop of the moments and the 
characteristic dip for p™* ~ 1 GeV indicate, however, that perturbative calculations may 
be relevant for hadrons down to a scale of approximately 1 GeV. 

The Monte Carlo model calculations which include both the parton cascade and hadro- 
nisation, largely follow the trend of the data and, in particular, reproduce the minimum 
around p™^ = 1 GeV. Differences appear for p^'^^ < 1 GeV, with Herwig describing 
the data better than the models using string fragmentation. However, the dotted curves, 
which represent Pythia predictions with the inclusion of Bose-Einstein correlations, are 
in very good agreement with the measurements. 

To study the influence of resonance decays, event samples were generated wherein 
resonances were not allowed to decay. The enhancement for p™^ < 1 GeV remained almost 
unaffected in the case of the Pythia model. For the Herwig model, the suppression of 
the decays led to an increase of the factorial moments below p™* ~ 0.5 GeV, an effect 
also observed in e''"p studies [12]. This increase ranges from a few percent for g = 2 to 
about 20% for g = 5. 

Fig. 121 shows spherically cut factorial moments Fg, g = 2 to 5, as a function of p™*. 
For large p™* values the moments change very little owing to the kinematically limited 
number of particles per event at high momenta. However, for smaller p*^"* the moments 
increase rapidly and show no tendency to level off for p"^"* — > GeV, contrary to theoretical 
predictions for partons. The Monte Carlo models describe the data rather well down to 
pcut ^ 2 — 3 GeV. In that region they are, in fact, very similar to parton-level (and hadron- 
level) predictions from Ariadne shown in [10]. For smaller values, and in sharp contrast 
to the data, the Monte Carlo curves flatten off and remain approximately constant, a 
feature also observed in [10]. The influence of Bose-Einstein correlations, as implemented 
in Pythia, is sizeable in that region but insufficient to reproduce the data. Indeed, 
the LEP measurements [33] have shown that BEG effect is expected to be larger in the 
longitudinal direction than that in the transverse direction of the jet, though this is not 
implemented in Pythia. 

4.2 Factorial cumulants 

The large statistics available in this analysis allows the study of the factorial cumulants, 
i^g(p™*) and Kg{p^^^), defined in (j2)). These are a direct measure of the genuine correla- 
tions among hadrons and thus present the information in a way which can more readily 
be interpreted. Fig. El shows the cylindrically cut and spherically cut cumulants of order 
g = 3,4.4 

The cumulant K^lp^^^^) has a similar p^^^^-dependence to the factorial moments: it is 
positive at large p™*, decreases towards a minimum, close to zero around 0.6 — 0.7 GeV 
and rises rapidly as p^^^ is further decreased. Interestingly, the minimum occurs at a 

"^Note that K2 is, by definition, equal to F2 — 1 and therefore not shown here. 



smaller p™* value than that for K2{p^"^^) (or F2 — 1, Fig. P) and for the higher-order 
factorial moments. Four-particle correlations, measured by K^lp™^), are compatible with 
zero, within errors. From the rapid increase of i^2(PT"*) cind K^lp^^) as p™* -^ GeV, 
we can conclude that the strong rise of the factorial moments is predominantly due to 
genuine particle correlations. 

Figure 01 shows that the cumulants in spherically cut phase space are significantly 
larger than the Kq[p^^) cumulants for p^^^ smaller than a few GeV. They continue to 
increase with decreasing p'^'^*, a behaviour reflected in the corresponding factorial moments 
in Fig. 

The various Monte Carlo models, shown in Fig. IHl agree qualitatively with the mea- 
surements but differ in detail for all p™* and p^'^^ < 2 GeV values. The largest deviations 
occur for the cumulants i^g(p™*) below p™^ ~ 2 GeV where the models start to level 
off, whereas the measured cumulants continue to increase as p™^ -^ GeV. Suppression 
of resonance decays in Herwig was found to increase the cumulants Kg{p^^) for p^"* 
below 1 GeV. In contrast, the spherically-cut cumulants predicted by Herwig remain 
essentially unchanged. 

The results presented above might potentially be biased due to correlations induced by 
Dalitz decays (vr^ —>■ 7e"'"e~) and fake pairs. The former would manifest themselves as a 
narrow peak in the invariant mass distribution of unlike-sign particle pairs near threshold. 
No such enhancement was found in the data. The check against fake pairs showed a peak 
at very small invariant masses which was found not to influence the results. 

We have also repeated the analysis for multiplets composed of like-sign particles. 
Although the corresponding moments differ in magnitude, their dependence on p™^ and 
pcut ^^^^ shown) follows closely that of all-charge particle moments. We may therefore 
conclude that the dip-structure observed in Fig. [T] and Fig. |21 is stable against changes in 
the charge composition of the multiplets. 

4.3 Comparison with deep- inelastic scattering 

The ZEUS measurements reported in [12] were carried out in the current region of the 
Breit frame of reference which is traditionally considered to be the equivalent of a single 
event-hemisphere in e"'"e~ annihilation [13]. A sample of e"'"p interactions was used with 
an average four-momentum transfer squared of {Q"^) ~ 2070 GeV^. This corresponds to 
an equivalent e"'"e~ c.m.s. energy, y/s, of 44 GeV. 

The factorial moments studied in [12] show many of the characteristics also reported 
here. However, an interesting difference is observed in the p^^^ dependence of the factorial 
moments. The distinctive minimum seen in the OPAL data for p™* ^ 1 GeV, which could 
signal the borderline between perturbative and non-perturbative dynamics, is absent in 
the ZEUS measurements. The latter remain constant down to p^^^ = 1 GeV below 
which value they increase rapidly. The authors interpret their measurements as the first 
indication that perturbative QCD fails on a qualitative level to describe the hadronic 
multiplicities and that hadronisation causes a violation of the LPHD hypothesis for many- 
particle inclusive observables. 
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To try and understand the differences between the e"'"e~ and e"'"p results, we have 
carried out further studies to mimic the ZEUS experimental conditions. These differ 
significantly from those in e"'"e~ annihilation: (i) the selected current region of the Breit 
frame excludes a large part of the central rapidity region in the 7*-proton rest-frame, the 
equivalent of the c.m.s frame in e"'"e~, which is included in our analysis; (ii) the rate of 
hard jet production is lower in e"'"p collisions at HERA energies than in e"'"e~ at LEP. 

To study the influence of the central rapidity region and the hard jet production rate, 
we have repeated the analysis for rapidity intervals y > i/o with yo '^ 1 and for two-jet 
events, the latter selected through a cut on the thrust value. 

As an example. Fig. |3] shows the moments Fq(p™*) for yo = 1.5. Omitting the central 
region in the OPAL data clearly has a strong effect on the magnitude and p™*-dependence 
of the moments. In particular, the minimum around p™* = 1 GeV, seen in the full sample 
{y > 0) is absent when only hadrons with y > 1.5 are selected. The ZEUS and OPAL 
data {y > 1.5) overlap for q = 2 in the full p™^ region; higher order moments are very 
similar below p™* of about 1 GeV but differ in magnitude at higher p™*. 

In Figure |3] we also show results for two-jet events. These events were selected by 
requiring the thrust value of an event to be larger than 0.96. The sample corresponds 
to about 30% of all multi-hadronic events. An alternative two-jet selection based on the 
Durham jet-finder [26] led to very small differences which are included in the systematic 
uncertainties. 

Selecting two-jet events reduces significantly the values of the factorial moments (and 
cumulants, not shown) for all p™* in comparison with those in the inclusive sample. The 
OPAL measurements are lower than the ZEUS data but, as expected, show qualitatively 
the same behaviour: constant down to p™* ~ 1 GeV followed by a strong rise as p™* -^ 
GeV, with no evidence for a distinct minimum near 1 GeV. 

The OPAL results on rapidity-restricted factorial moments suggest that the often 
assumed equivalence of a single event-hemisphere in e+e~ annihilation with the current 
region in the Breit frame for deep-inelastic e"'"p interactions has to be treated with caution. 
In theoretical predictions for QCD cascades, which apply to the whole jet and moreover 
focus on soft gluon emissions, characteristic signatures of soft particle emission, such 
as colour coherence, are seen only if particles with small rapidities are also included. 
Likewise, selecting two-jet events may introduce a bias which can mask the effect under 
study. Our findings suggest an explanation for the differences obtained in the ZEUS 
analysis [9] of angular correlations compared to the LEP results [8]. 

5 Summary and conclusions 

Analytical perturbative QCD calculations, in agreement with parton-level Monte Carlo 
calculations, show that gluons produced in a jet become uncorrelated when the gluon 
transverse momentum relative to the jet axis, px, is restricted to small values. The 
approach to a Poisson regime is a direct consequence of colour coherence, or angular 
ordering of gluon emissions in the QCD cascade and is expected to hold also for soft 



hadrons if Local Parton-Hadron Duality (LPHD) is valid. 

In this paper, the predicted QCD colour coherence effect has been tested in e"'"e~ 
annihilation at the Z-resonance using factorial moments and factorial cumulants of the 
multiplicity distribution of hadrons with restricted transverse momenta pt < Pt"* or 
restricted absolute momenta p = |p| < p^^^. The analysis is based on a data sample of 
about four million events recorded with the OPAL detector at LEP. 

For cylindrically cut phase space, the factorial and cumulant moments are predicted by 
analytical QCD calculations to reach limiting values close to unity and zero, respectively, 
for small p™*. Likewise, in spherically cut phase space, the factorial and cumulant mo- 
ments should saturate at values well above unity and zero, respectively, for p'^'^* — * GeV. 
These expectations are not borne out by the measurements: for p|."* < 1 GeV or for values 
Q^ pcnt i^giQ-^ a few GeV, the moments and cumulants rise strongly with decreasing p™* 
or p'^^^, QCD based Monte Carlo models which include hadronisation, reproduce well the 
change in correlation pattern but serious discrepancies remain in the very low momentum 
region. However, Bose-Einstein correlations, as implemented in Pythia, significantly 
improve the agreement between data and Monte Carlo predictions. 

Interestingly, in the region of large to intermediate p™*, a minimum around 1 GeV is 
observed in the p^^ dependence of the factorial moments, with a corresponding minimum, 
close to zero, at p^^ ~ 0.6 GeV for the cumulant K^, as expected from colour coherence. 
One may interpret the intermediate p™*-range 0.6 — LO GeV as a borderline between a 
regime of perturbative dynamics where the LPHD hypothesis is justified, and a regime 
dominated by strong confinement forces which leads to violation of LPHD for many- 
particle inclusive observables. 

In a similar analysis in deep inelastic e"''p scattering at HERA, no evidence was found 
for a Poisson-like regime in cylindrically-cut phase space. The results presented here show 
that the characteristic decrease towards a minimum around p^"* of 1 GeV disappears if 
hadrons produced in the central region of rapidity (y < 1) are excluded, or if the analysis 
is restricted to two-jet events. The former suggests that, for soft particle production, 
the often assumed equivalence of a single event-hemisphere in e"'"e~ annihilation with the 
current region in the Breit frame of a deep inelastic e"'"p collision may be misleading. 
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Figure 1: Factorial moments of charged particles with transverse momenta j9t < p™^ as 
a function of p^^ compared to different Monte Carlo predictions. The error bars, shown 
where larger than the marker size, represent the statistical and systematic uncertainties 
added in quadrature. The dotted line shows Monte Carlo predictions from Pythia with 
Bose-Einstein correlations included (see text). The errors of the Monte Carlo predictions 
are comparable to those of the data. 
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Figure 2: Factorial moments of charged particles with absolute momenta p = |p| < p^'^^ as 
a function of the the momentum cut, p'^"* compared to different Monte Carlo predictions. 
The error bars, shown where larger than the marker size, represent the statistical and 
systematic uncertainties added in quadrature. The dotted line shows Monte Carlo predic- 
tions from Pythia with Bose-Einstein correlations included (see text). The errors of the 
Monte Carlo predictions are comparable to those of the data. 
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Figure 3: K^ and K^^ cumulants of charged particles with transverse momenta px < Pt"* 
as a function of the momentum cut p^^ (left panel) and with absolute momenta p < p™* 
as a function of the the momentum cut p'^"* (right panel) compared to different Monte 
Carlo predictions. The error bars, shown where larger than the marker size, represent the 
statistical and systematic uncertainties added in quadrature. The errors are correlated bin- 
to-bin. The dotted line shows Monte Carlo predictions from Pythia with Bose-Einstein 
correlations included (see text). The errors of the Monte Carlo predictions are comparable 
to those of the data. 
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Figure 4: Factorial moments of charged particles with absolute momenta pt < Pt"* '^^ '^ 
function of p^^ , compared to those of 2- jet events, in the rapidity window y > 1.5 and 
data from ZEUS [12]. The error bars, shown where larger than the marker size, represent 
the statistical and systematic uncertainties added in quadrature. 
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